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Abstract 

We present the results of a high statistics lattice study of the gluon 
propagator, in the Landau gauge, at (3 = 6.0. As suggested by previ- 
ous studies, we find that, in momentum space, the propagator is well 

described by the expression G(k 2 ) = M 2 + Z ■ k 2 [k 2 / A 2 ) 11 

By comparing G(k 2 ) on different volumes, we obtain a precise de- 
termination of the exponent r/ = 0.532(12), and verify that M 2 does 
not vanish in the infinite volume limit. The behaviour of r\ and M 2 
in the continuum limit is not known, and can only be studied by in- 
creasing the value of (3. 
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Introduction 



One of the most interesting, yet unresolved, issues concerning Quantum Cro- 
modynamics, as a theory of strong interactions, is the comprehension of 
the transition from the fundamental fields, appearing in the Lagrangian, to 
the physical particles, namely the hadrons. While color confinement and 
hadronization are believed to occur, no consistent derivation of this mecha- 
nism from first principles is yet available, although much effort has been put 
in this direction. 

By studying the propagation and the interaction of the fundamental 
fields, appearing in the QCD Lagrangian, one attempts to gain some in- 
sight on the mechanism of confinement and hadronization, which are gen- 
uine non-perturbative effects. In this paper, we present the results of a non- 
perturbative study of the gluon propagator, in the Landau gauge, obtained 
from a high statistics lattice simulation. 

The results of the present study, on a volume larger than in refs. Jl], @], 
confirm that the gluon propagator is well described by the functional form 
proposed in the references above: 



We obtain a more precise determination of the exponent, finding t] = 0.532(12) 
and show that the gluon propagator at zero momentum stays finite, as the 
infrared cut-off is removed. 

In section 1, we review what is known from perturbation theory, and some 
results obtained by perturbative analyses. Particular attention is devoted to 
other lattice simulations, whose results can be compared directly with ours. 
In section 2, the definition of the correlation functions will be given, together 
with a description of the numerical implementation. Section 3 is dedicated to 
the discussion of systematic effects and lattice artefacts. We found important 
to discuss these aspects before giving the results, which are largely influenced 
by our control of systematic errors. Finally, in section 4, we present the new 
results and compare them with those of our previous study 0, trying to 
interpret common features and differences. 
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1 The gluon propagator in the literature 



As is the case for gauge-dependent quantities, the gluon propagator can only 
be defined once that the gauge has been fixed. In the following, we will 
solely be concerned with the Landau gauge, which corresponds to fulfill the 
condition 

d lt Afa) = 0, a = l,...,8, (2) 
where A a ^(x) are the gauge fields, and a the color index. 
The gluon propagator in this gauge has the form 

Dftv(k) = -i[ d*x(\T (A,(x)A u (0)) |)e fa = G(k 2 ) (g, u - ^] , (3) 



k 2 t 

where G(k 2 ) = —1/k 2 at tree level. When leading order perturbative correc- 
tions are taken into account, G(k 2 ) takes the form 

G{k ' ^ = ~k 2 (l + U(k 2 )) ~ ~(k 2 y+n> ( 4 ) 

with rf ert = — 7g(7 2 /47t and 7g = —3/13. Therefore, at values of the coupling 
constant corresponding to (3 = 6.0, g 2 ~ 1, these corrections correspond to a 
minor change, at most 5%, to the tree level expression^ 

Perturbative corrections to G(k 2 ) cannot therefore account for the higher 
power of k 2 which is necessary to give a confining potential, in the one gluon- 
exchange approximation. 

An analytic non-perturbative approach to the problem is indeed very 
difficult. The most elegant and neatest way to attack the problem is per- 
haps to solve the Dyson-Schwinger equations w.r.t. the renormalized gluon 
two-point function D^k). The infinite set of differential equations should 
satisfy the Slavnov- Taylor identities in order to preserve gauge invariance. 
Being the equations non linear, different solutions are allowed 0-0, and 
their interpretation, already difficult, is made even more uncertain by the 
approximations that one has to make to render the problem feasible. 

1 To obtain this extimate we have not considered the running of the coupling constant. 
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In ref. || a procedure for deriving the gluon propagator in the Minimal 
Landau Gauge, is presented. Subsequently, Zwanziger confirmed the so- 
lution of ref. by repeating the calculation on the lattice. In ref. [|IJ, we 
have shown that this expression fails to reproduce the lattice behaviour. 

In ref. || , the first exploratory study on the lattice was presented. They 
found that the effective energy is an increasing function of the time (see below 
for definition), that the gluon propagator violates the spectral decomposition 
and, in particular, cannot have the form (|3]) with G(k 2 ) ~ 1/ (k 2 +m 2 ). Subse- 



quently, the authors of ref. || (TOJ] tried to account for such a behaviour either 
by including the possible contribution of unphysical particles, or considering 
the solution proposed in refs. || [7|. 

Lately, to interpret the results of the lattice simulations, a new functional 
form has been proposed [ffl, 0]: 



G{k 2 ) 



M 2 + Z-(k 2 )l^ 2 



fc2\T|- 1 



(5) 



This expression tries to mimic the effect of an anomalous dimension for the 
gauge field, and keeps G{k 2 ) finite. The latter is a necessary condition, due 
to the fact that the lattice possesses an intrinsic infrared cut-off. 

According to |TJ r) — 0.4 0.6, whereas the authors of |2j quote r\ = 
0.35. The value of the dimensionful parameter M 2 is more controversial, and 
difficult to determine because it is very small, of the order of the infrared 
cut-off. It is not clear whether M 2 will vanish in the infinite volume limit. 

Since there is no a priori theoretical prediction of the functional form (|5|), 
several further investigations are in order. Namely, 

1. reduce the ultraviolet cut-off by increasing the value of (3, and try to 
determine M 2 and r\ in the continuum limit; 

2. verify if (|5|) contains any gauge-independent information. This can be 
investigated by studying the gluon propagator in different gauges. 

In this paper, we address a different issue, that is the dependence of 
the parameters of equation ([5]) on the infrared cut-off. Anticipating some of 
the results, we have indications that the dimensionful parameter M 2 in not 
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merely a finite volume artefact, and it would keep a non-zero value in the 
infinite volume limit. Its fate in the continuum limit, i.e. as (3 — > oo, is not 
known. The value of the anomalous dimension r], on the other hand, does 
not seem to be affected by the presence of the infrared regulator. A major 
question is the interpretation of a finite value of the gluon propagator, at 
zero momentum. 

2 The gluon propagator on the lattice 

In this section, we establish the notation and the main formulae, to be used 
in the following. We also give a brief description of the algorithms, i.e. the 
Monte Carlo procedure, used to generate statistically independent configu- 
rations, and the method which has been used to fix the Landau gauge. We 
conclude by reporting some information about the accuracy achieved by our 
gauge-fixing algorithm. 

2.1 Notations 

The gluon propagator is studied on the lattice by looking at the euclidean 
two-point correlation functions: 



where /i, v = 1,...,4 and the trace is intended over color indices. The 
corresponding expression in momentum space is 



Among the possible definitions of the gluon field on the lattice, we choose 



In a previous publication [Q , we have discussed in detail the motivations for 
such choice. 




(6) 






(8) 
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In order to study the gluon propagator in momentum space, D^ u (k), we 
take the 4-dimensional Fourier transform of the correlators Taking into 
account translation invariance, is given by 

Dfjcik) = J] Tr(A^x)A u (0))e lkx . (9) 

x,t 

On the lattice, k can only take discrete values, k^ = (2ir/ L^ajn, where a is 
the lattice spacing, and n = 0, 1, . . . , L M — 1. 

Moreover: 

• we always take full advantage of space-time translational symmetry by 
averaging the previous expression over different points on the euclidean 
lattice; 

• from r.h.s. of @ the scalar function G(k 2 ) is extracted, by taking suit- 
able combinations of the spin indices //, v accordingly to (|7|). 

• in order to check the validity of continuum dispersion relation on the 
lattice, the following four classes of correlators have been considered: 

1. Transverse correlations 

D u (k) = J2^(M^M0))e ik ^e ikot = G(A; 2 )4, (10) 

x,t 

where i, j = 1,...,3 and the spatial momentum kj is non-zero 
only along axis orthogonal to the spin direction i. 

2. Longitudinal correlations 



Djj(k) = Y / Tr(A J (x)A J (0))e^e^ t = G{k 2 ) ^ - , 

(11) 

where j = 1,...,3 and the spatial momentum has a non-zero 
component along the direction of the spin j. 

3. Non-diagonal correlations 

Dij(k) = '£Tr(A i (x)A j (0))e il3 e ikot = -G(k 2 )^, (12) 
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where i,j = 1, . . . , 3 and % 7^ j. These correlations do not vanish 
only when the spatial momentum is non-zero along both % and j 
directions. 

4. Temporal correlations 

D 00 (k) = £Tr(A (x)A (0))e l H^ = G{k 2 ) U m - M) . (13) 

x,t \ / 

Further investigations of the behaviour of the gluon propagator can be 
carried on by looking at some correlations, which are functions of the spatial 
momentum and euclidean time. For simplicity we maintain the same notation 
given before, namely, we will consider: 

1. Transverse correlations 

D«{t, k) = J2^(Mx)M0))e ik ^, (14) 

X 

2. Temporal correlations 

D 00 (t,k) =£Tr<A o (z)A(0)>e^. (15) 

X 

All the previous expressions have been computed for the following val- 
ues of spatial momenta: k = (0, 0, 0), (27r/24a, 0, 0), (47r/24a, 0, 0), and k = 
(27r/24a, 27r/24a, 0). Moreover, permutations of the three spatial axes have 
been considered. The Fourier transform, with respect to the time coordi- 
nate has been computed allowing k$ to take all the possible values. Namely, 
k = (2n/A8a)n, with n = 0, . . . , 47 . 

2.2 Lattice setup and gauge-fixing 

We summarize the parameters of the present simulation in tab. [I]. The only 
differences, with respect to ref. Q, are the lattice size and the number of 
configurations. All the others parameters of the simulation have been left 
unchanged, with the purpose of studying the dependence of expression (H) on 
each of them. For the generation of the pure gauge fields we adopted a Hybrid 
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P 


# confs. 


Volume 




present work 


6.0 


500 


24 3 x 48 


< 10~ 6 


ref. § 


6.0 


1000 


16 3 x 32 


< io- 6 



Table 1: Summary of the parameters of the present simulation and of those 
of the simulation of ref. [1] 

Monte Carlo (HMC) algorithm in its local version and with a checkerboard 
update, as described in |TT|. This version of the HMC allows very short 



autocorrelation times among pure gauge configurations and is very effective 
when high statistic measures are needed. In our simulation we performed an 
initial termalisation of 1024 iterations and then we extracted configurations 
every 128 iterations. 

A further comment on the size of the statistical sample: although the 
number of configurations generated during the present run, is a half the old 
one, we have obtained more accurate results. This is due to the fact that 
the larger number of lattice points (the ratio of the two volumes is about 8) 
improves the stability of the signal. 

In presence of gauge-dependent quantities, the accuracy of the gauge- 
fixing algorithm is a crucial information. Uncertainties in the gauge-fixing 
result in residual fluctuations, thus causing extra noise and/or systematic 
alterations. For this reason, the accuracy can only be established by per- 
forming tests on the gauge-dependent operators under study. 

As demonstrated in ref. [l[] , the residual fluctuation left-over after gauge- 
fixing 

(^(a^Latt < 10" 6 , (16) 

are absolutely negligible, with respect to the statistical errors. Since we 
consider this to be a crucial point, we have repeated the same analysis as in 
[ l]] on the new sample. We have defined logaritmic (symmetric) derivatives, 
with respect to t, of transverse and temporal correlators, at zero spatial 
momentum. As in ref. PJ, we find that the temporal correlator is constant 
and its derivative always compatible with zero, within errors. This qualitative 
behaviour can be compared with a similar figure in ref. M. 
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t D 00 (t,0) 


t D 0O (t,0) 


t Ax)(*,0) 


1 56552 

2 56552 

3 56550 

4 56549 

5 56548 


6 56547 

7 56549 

8 56548 

9 56547 
10 56548 


11 56545 

12 56547 

13 56542 

14 56543 

15 56543 



Table 2: Doo(t, 0), measured on one configuration, as a function of the time 
separation. The gluon field A (t,k) is constant in time, as required by the 
gauge condition. 

A different check can be performed by noting that the condition (0) im- 
plies that the three dimensional Fourier transform of the gauge field Ao(t, k = 
0) is time independent. This time independence should hold on each single 
configuration. This has been checked by looking at the temporal correlation 
D o(t, 0) of some arbitrarily chosen configurations. In table |2], we report the 
value of D 00 (t,0), at different time separations, on a single configuration. 
Doo(t,0) is constant at the level of 0.008%, which is in fact consistent with 
the precision dl6|) . Results in table |]are to be compared with figure 1 of ref. 
[0], where the gluon field, plotted versus t after gauge-fixing, shows a strong 
time dependence. 

The previous discussion shows that the present gauge-fixing procedure is 
the most effective one, among those implemented in the literature. 

3 Results for the gluon propagator 

In this section, we outline the unusual behaviour of the effective energy. 
Then, we consider the dependence of G(k 2 ) on the momentum. Assuming 
the validity of the functional form (|5|), we first discuss in length systematic 
effects, and give our final results. 
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3.1 The behaviour of the effective energy 



Using spectral decomposition and translation invariance, the correlation func- 
tions can be rewritten as 

fl (UVE«e- (17) 

|<) ^ 

where the sum is intended over all the states with the quantum numbers 



of operator A(x) and Ei = ymf + \k\ 2 is the energy of the state i. The 
normalization Mi is positive for physical states. From eq. (|P7D, the effective 
energy 

u; eS (t,k)= log J^'fl - (18) 
D(t + 1, k) 

should be a decreasing function of the time, for any value of the momentum 
k. 

In figure [I], we report 0J e fi{t, k) computed from the transverse correlators 
(p!4|). uo e s {t, k) is clearly increasing with time, for all the momentum combi- 
nations considered. Some authors || tried to justify this result by including 
( |T7D the contribution of unphysical particles, whose normalization is nega- 
tive. The function apt to fit the lattice points, at large time separations, 
would then be: 

F(t) = ^e- mit + -^e" m2i ; Z 1 < 0, Z 2 > 0; m x > m 2 . (19) 
Ami 2rri2 



In our previous paper |IJ, we showed that fll^D fails to reproduce the be- 
haviour of D^it, k). The apparent agreement of flT9"| ) with D(t,k = 0) in ref. 
[0], was probably faked by the large statistical errors. 



3.2 Systematic effects 

The systematic effects discussed in the following are those due to the infrared 
and ultraviolet cut-offs, introduced when the theory is regularized on a finite 
lattice. These effects limit the range of momenta, where reliable results can 
be obtained. A way to detect such effects is to look for discrepancies between 
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continuum dispersion relations and their lattice realization. We note that, 
contrary to what happens for other quantities, the quenched approximation 
is not going to produce systematic changes, being the pure gauge sector a 
self-consistent part of the theory. One then expects these results to stay 
valid per se, as a description of the pure gauge behaviour, although they 
might differ in the full theory. Another source of systematic uncertainty is 
the value of lattice spacing, which is found to vary substantially depending 
on the physical observable used to fix it. We will only be concerned about 
this point at the end of this section, when we will express, in physical unit, 
the propagator at zero momentum. 

We discuss the effects of the ultraviolet cut-off, first. It is known that 
the presence of a finite ultraviolet cut-off produces unwanted effects, when 
momenta are of the order of the cut-off itself. On the lattice, this means that 
results, corresponding at momenta such that k 2 a 2 > 1, contain discretization 
contributions. A demonstration of the above statement, is given in figure 
|^A, where G(k 2 ) is plotted as a function of k 2 in dimensionless units. In 
correspondence of the same value of k 2 , G(k 2 ) can have different values, due 
to the breaking of the Lorentz symmetry. The different results correspond to 
G(k 2 ), as extracted from the four classes of correlators (|10|)-(|T3|). They agree 
within errors until k 2 a 2 ~ 1 where they start clustering into two separated 
curves. This behaviour, at larger momenta, is caused by the fact that, in 
the Fourier transforms (|10|)-(|l"3|), we only increase k Q , to the larger values 
allowed on the lattice, while k iy i = 1,2,3 are always kept small. Points lying 
on the higher curve are those generated through spatial correlators (|iT^)-(|T2|), 
while, on the lower curve, there are the points computed from the temporal 
correlators (|I3|), which are directly proportional to /c 2 and so more affected 
by discretization effects, at large k . 

To verify this hypothesis, we have tried a kinematical correction to the 
momentum, substituting 

k, — /C M = 2sin(^), 

k 2 — /C 2 = E(2sin(|i)) 2 (20) 

in the dispersion relation (|7]). Such a correction is suggested by the tree 
level expression of the bosonic propagator on the lattice. In fig. QB, G(k 2 ), 
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obtained by performing the substitutions (p0[), is shown. One sees that the 
two curves of fig. have collapsed into a single one. On a logarithmic 
scale, the curve, obtained by using ( f2"0| ) appears as a straight line, at large 
k 2 , whose slope is the exponent 1 + rj of equation ([5]). Around JC 2 a 2 = 1 
the curve changes slope and rj is slightly reduced. We cannot tell whether 
this behaviour is due to residual cut-off effects, or to the fact that, as the 
momentum increases, rj is going towards its perturbative value. This issue 
can only be decided by decreasing the lattice spacing, i.e. by going to larger 
values of (3, and increasing the range of k 2 , where lattice artefacts are small 
and momenta are large in physical units. 

The investigation of the effects of the infrared cut-off is delicate and is 
directly related to the determination of the value of M 2 in ( At very 
low momenta, different determinations of G(k 2 ) corresponding to the same 
value of k 2 are inconsistent (see fig. 3). This effect is due to the fact that 
our lattice has different length, along the spatial and temporal axes. In fig. 
3, we indicate, with an arrow, the value of k 2 where the differencies became 
negligible. 

We deduce from the previous analysis that one should only take into 
consideration the region of of intermediate momenta. More quantitative 
indications in this sense are provided the fitting procedure, see below. 



3.3 Fitting procedure 

In the previous section, it has been shown that infrared and ultraviolet effects 
are reasonably small, only in the intermediate region of momenta. On large 
volumes, the region of momenta "free" from cut-off effects contains many 
points, which can be used to check the stability of the fits. The fitting 
function is the dimensionless version of (||) 

G( * 2) = Ml + Z L (,^ (21) 

which depends on three free parameters. We fit first G(k 2 ) to expression 
(pip, on an interval of 6 points, which we move back and forth in the selected 
region of momenta. The fit is then repeated using the same procedure with 
8 and 12 points, respectively. In the central region of k 2 , the results are 
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perfectly compatible on different fc 2 -ranges and their stability improves, by 
increasing the number of points. 

In figure we plot the best 12-points fits for 7] and M 2 , together with 
the corresponding value of xidof> as a function of the minimum momentum 
considered in the fit (k^). As anticipated in the previous section, for low 
and high momenta the fits are rather poor, giving unstable values for the 
parameters and a large Xndof- At momenta k^^a 2 ~ 0.1 the values of Xndof 
lower to 1 -r 1.5 and the estimates for 77 and M 2 become very stable. Our 
final determination obtained by averaging over the four points indicated in 
the figure, is 

Ml = 4.46(9) x 1(T 3 
Zl = 0.102(1) 
T] = 0.532(12) 
(xLof) = i-os 

One can compare the present results with those of ref. Jj]]. We have re- 
analyzed the old result, by fitting G\{t) to expression (^1])Q. On a volume 
V = 16 3 x 32 we find: 

Ml = 2.8(1) x 10~ 3 
Z L = 9.01(4) xlO- 2 

r] = 0.56(6) 1 ' 

Xndof = 1-5 



4 Discussion of the results 



Apart from the value of M 2 (see next section), the agreement between the two 
sets of results is reasonably good, and the quality of the fits is satisfactory. 
By increasing the volume, we have obtained a more accurate estimate of the 
exponent 77. 

In our previous paper, we have quoted three different sets of results, 
obtained by studying different correlators, as functions of k 2 or k, t respec- 
tively. We interpreted the differences as systematic effects. We are now able 
to check this interpretation. On a smaller volume, we had a limited choice 

2 In is necessary to change the overall normalization of G\(t), in order to make it 
independent on the lattice volume. 
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of momenta, before running into problems when k 2 a 2 ~ 1 or infrared effects 
become important. As figure |] shows, some of the points we had considered 
in ref. [TJ, in the fit of the propagator in momentum space, (see eq. (19) in 
ref. HI), correspond to incompatible estimates of G(k 2 ). For this reason, in 
ref. JTJ , the fit performed in momentum space, was giving results uncompat- 
ible with those of eqs. (22) and (23), and it also disagrees with the present 
result fl22|) . 

In ref. [0], all the values of momenta, (including those corresponding to 
k 2 a 2 > 1) have been considered in the fit. Due to the flattening of the curve 
in the region of large momenta, the value of r] has been underestimated. 

Finally, we observe that the finite value of V has a very small effect in 
the determination of rj, provided that the range of k 2 is large enough. 



4.1 Comments on the value of M 2 

The main difference between (22) and ( ^3l) is the value on the parameter M 2 , 



which increases with the volume. This observation rules out the hypothesis 
that the non-zero value found for M 2 is merely due to finite volume effects. 
If this were the case, we would expect M 2 to scale roughly as 1/L 2 . 

Further arguments may be provided to reinforce the conclusion that 
M 2 ^ 0. In figure |3], values of G(k 2 ) measured on a 16 3 x 32 lattice are 
over- imposed to the present data. We notice that G(k 2 = 0), whose inverse 
gives an estimate of M 2 , is smaller on the largest volume. This is opposite 
to the expectation M 2 oc 1/L 2 . 

A second argument is provided by the comparison of the fit to the data. 
In figure ^, the best fitting curve is compared to the lattice numbers, showing 
a good agreement even outside the fitting range. On the right side, the ratio 
data to fit shows that the fit overshoots the data, at low momenta. That 
is, the behaviour of the propagator at the so called intermediate momenta 
is consistent with a value of M 2 larger than G(k 2 = 0) found on the lattice. 
This is a further proof that the effect of the infrared cut-off is to fake a value 
of M 2 lower that the real one. 

It is possible to try a first, very crude extrapolation of the value of M 2 
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to the infinite volume limit. Using 



M 2 {V) 



M 2 (V = oo) + cost 



1 



(24) 



y/V' 



we get 



M 2 (V 



oo) = 6.202(8) x 10 



-3 



(25) 



By using a 1 ~ 2 GeV, as it is determined by several lattice simulations 
at (5 = 6.0, this corresponds to 



It is fundamental to understand the behaviour of M 2 in the continuum, 
i.e. as (3 — > oo. Notice that, even if M 2 stays finite in this limit, it cannot be 
interpreted as a symptom of symmetry breaking, since, in (|^), M 2 does not 
represent a pole. An interpretation of this result, in connection with colour 
confinement is, at present, absent. 

5 Conclusions 

The study of the gluon propagator on different lattice volumes at fixed value 
of the coupling constant, confirms that G(k 2 ) is well described by the ex- 
pression ([5]). The value of the anomalous dimension rj has been confirmed 
to assume a value slightly exceeding 0.5. By increasing the volume from 
16 3 x 32 to 24 3 x 48, we have obtained a better determination of the anoma- 
lous dimension: rj = 0.532(12). We find that M 2 does not vanish as 1/L 2 . 
Assuming a -1 ~ 2 GeV, our best estimate is Mp hys ~ (160MeV) 2 . 

Further investigations to study the continuum limit and the gauge-dependence 
of the results presented above are certainly required. If, ultimately, M 2 
should stay finite, the comprehension of the mechanism which allows the 
gluon to develop a mass-scale in the infrared region would call for a theoret- 
ical understanding. 



Ml 



phys — 



(160 MeV) 2 ~ (Aqcd) 2 



(26) 
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u cS (t, k) 



t 

Figure 1: Effective energy for transverse correlators (20). The effective en- 
ergy is increasing with time for all the values of k. The four curves cor- 
respond to the following momenta: □ : k = 0, O : k = (2n/24, 0, 0), 
x: k = (2tt/24, 2tt/24, 0) and : k = (4vr/24, 0, 0) 



16 



G-\k 2 ) 

Figure A 



Figure B 



k 2 [a 2 ] IC 2 [a 2 } 

Figure 2: Behaviour of G^ 1 {k 2 ) at high momenta. Fig. 2A: we observe the 
data splitting into two curves at high momenta (k 2 a 2 ~ 1). Fig. 2B: the two 
curves are mapped into a single line when the kinematical correction (26) is 
applied. The logarithmic behaviour of the corrected curve is much closer to 
a straight line. 
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V = 16 3 x 32 

V = 24 3 x 48 



k 2 \a 2 



Figure 3: Behaviour of G(k 2 ) at low momenta. The arrow indicates the value 
of k 2 where different estimates of G(k 2 ) start to agree, within the errors. 
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k 2 [a 2 ] k 2 [a 2 } 

Figure 4: Best fit for the propagator in momentum space with the function 
(27). The figures show the fit (left) and the data to fit ratio, the dashed lines 
indicate the fitting range. The ratio shows that the fit overshoots the data 
for small momenta. 
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Figure 5: Values for M 2 and 77 versus the minimum momentum considered. 
The correspondent Xndof i s a l so shown. In this fit we considered 12 points. 
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